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William  D.  Davis,  Thomas  Cleary,  Michelle  Donnelly,  and  Samuel  Hellerman 

Abstract 

An  algorithm  for  predicting  smoke  or  carbon  monoxide  concentration  in  the  ceiling  jet  in  the 
presence  of  a smoke  layer  is  developed.  The  analysis  will  be  restricted  to  fires  that  produce 
turbulent  plumes  and  can  be  represented  by  axisymmetric  point  sources.  The  predictions  of  the 
algorithm  for  smoke  and  carbon  monoxide  are  compared  with  experiments.  Once  the  smoke  or 
carbon  monoxide  concentration  in  the  ceiling  jet  is  predicted,  the  activation  times  for  detectors 
can  be  calculated  using  a model  for  detector  activation. 

Introduction 

Predicting  smoke  detector  response  to  a growing  fire  requires  calculating  the  time  dependent 
evolution  of  the  smoke  concentration  in  the  ceiling  jet.  Typically,  the  temperature  rather  than 
the  smoke  concentration  is  used  to  predict  smoke  detector  response  due  to  the  availability  of 
correlations  that  give  ceiling  jet  temperature1  and  the  assumption  that  smoke  concentration  can 
be  related  to  ceiling  jet  temperature-’  . Using  temperature  to  predict  smoke  detector  activation 
ignores  differences  in  the  production  of  smoke  by  burning  materials  that  may  completely 
invalidate  a temperature/smoke  prediction  correlation. 

There  have  been  efforts  to  use  computational  fluid  dynamic  (CFD)  methods  to  calculate  the 
smoke  concentration  in  the  ceiling  jet  and  with  the  increased  computer  power  available  today, 
these  methods  are  becoming  practical4,5.  However,  there  is  still  a need  for  an  algebraic 
correlation  that  would  yield  smoke  concentration  in  the  ceiling  jet  and  not  require  substantial 
computer  power  to  obtain  the  solution.  Early  work  along  this  line  can  be  found  in  Alpert ' s 
paper  on  the  ceiling  jet  that  resulted  in  the  successful  unconfined  ceiling  jet  temperature  and 
velocity  correlations6  in  use  today.  Later,  Yamauchi7  extended  Alpert ' s work  to  calculate  the 
smoke  concentration  and  smoke  detector  activation  in  the  ceiling  jet  when  a hot  layer  was 
developing.  Yamauchi ' s method  required  the  solution  of  a set  of  differential  equations  in  order 
to  define  the  ceiling  jet  properties  as  well  as  a zone  model  to  define  the  depth  and  temperature 
of  the  hot  layer. 

In  this  paper,  an  algebraic  correlation  for  smoke  or  carbon  monoxide  concentration  in  the 
ceiling  jet  in  the  presence  of  a smoke  layer  will  be  developed.  The  analysis  will  be  restricted  to 
fires  that  produce  turbulent  plumes  and  can  be  represented  by  axisymmetric  point  sources. 

Once  the  smoke  or  carbon  monoxide  concentration  in  the  ceiling  jet  is  predicted,  the  activation 
times  for  detectors  can  be  calculated  using  a model  for  detector  activation. 
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Theory 


The  development  of  a ceiling  jet  algorithm  requires  the  solution  of  three  separate  problems. 

The  first  problem  is  to  model  the  smoke  plume  as  it  rises  from  the  fire  to  the  ceiling.  The 
second  problem  is  to  model  the  turning  point  as  the  plume  flow  impinges  on  the  ceiling  and 
turns  to  form  the  ceiling  jet.  The  third  problem  is  to  model  the  ceiling  jet  flow  as  it  slows  and 
deepens  in  its  movement  across  the  ceiling. 

The  Plume  Region 

In  order  to  develop  an  algebraic  correlation  for  the  smoke  or  gas  concentration  in  the  plume,  the 
following  assumptions  must  be  made  in  order  to  simplify  the  equations.  The  fire  will  be 
represented  by  a point  source  and  assumed  to  be  axisymmetric.  The  zone  model  approximation 
of  homogeneous  temperature  and  particle  densities  in  each  layer  is  assumed.  Gaussian  shapes 
in  the  radial  direction  will  represent  the  velocity,  temperature,  and  smoke  profiles  in  the  plume. 
All  air  entrained  into  the  plume  below  the  smoke  layer  interface  will  be  considered  to  be  smoke 
or  gas  free. 

The  mass  flux  of  smoke  or  gas  in  a radially  symmetric  plume  can  be  written  as 


ms(z)  = \Csp{r,z)uz(r,z)2nrdr 
o 
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where  Csp(r,z)  is  the  scalar  representing  the  mass  concentration  of  smoke  particles  or  a gaseous 
combustion  product  such  as  CO  in  the  plume,  uz(r,z)  is  the  plume  velocity,  r is  the  radial 
distance  from  the  plume  centerline  and  z is  the  height  above  the  fire  source.  The  assumed 
Gaussian  profiles  for  the  scalar  mass  concentration  in  the  plume  and  the  plume  velocity  are, 


Csp{r,z)  = Csp0(z)e 
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where  Cspo  is  the  scalar  mass  concentration  for  the  plume  centerline,  uzm  is  the  plume  centerline 
velocity.  Act  is  the  1/e  width  of  the  plume  scalar  profile  and  a is  the  1/e  width  of  the  velocity 
profile.  It  has  been  assumed  that  the  smoke  or  gaseous  combustion  product  profile  in  the 
plume  scales  with  the  temperature  profile  in  the  plume2.  Integrating  equation  1 and  solving  for 
the  maximum  scalar  mass  concentration  at  the  plume  centerline  gives 
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The  plume  correlations  developed  by  Heskestad8  for  plume  radius  and  plume  centerline  velocity 

b and  uzm,  where  b is  the  plume  radius  where  the  temperature  has  dropped  to  0.5  the  plume 
centerline  temperature  T(z)  and  zo  is  the  location  of  the  virtual  point  source  with  respect  to  the 
fire  surface,  are  given  by 


Hz)  = C, 
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The  plume  centerline  scalar  mass  concentration  then  becomes 
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where  the  distance  z-zo  is  the  distance  above  the  virtual  point  source.  The  constant  C]  = 0.12, 
and  the  plume  radius  for  the  temperature  to  drop  to  1/e  of  its  centerline  value,  assuming  a 
Gaussian  distribution,  g-  1.201  b are  experimentally  determined. 

The  mass  flux  of  the  scalar  produced  by  the  fire  may  be  calculated  using  the  data  provided  by 
Tewarson9 


• v Q 

= 1 V 
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where  hc  is  the  heat  of  combustion,  Y/  is  the  scalar  yield  fraction  and  Q is  the  heat  release  rate. 
The  plume  centerline  scalar  concentration  is  given  by 


DQ *2/3 
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and  Z-z-Zo • The  temperature  ratio,  T(z)/T can  be  evaluated  using  Heskestad ' s correlations 
for  plume  centerline  temperature  excess8,  7fz)  - 7^,  and  location  of  the  virtual  point  source. 
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z0  = - 1.02  d + 0.083(22'5 
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where  d is  the  fire  diameter.  Hence, 


l + KQ 
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+ Smoke  yield  fraction  (in  grams  of  smoke  produced  per  gram  of  fuel  burned)  is  tabulated  in  many 
literature  sources.  Smoke  yield  fraction  can  also  be  estimated  by  dividing  specific  extinction  area  (from  the  cone 
calorimeter  ASTM  1354)  by  8.71  x 103  m2/kg10. 
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where  K=9.1(l-%rf/3 . 


The  plume  centerline  scalar  concentration  becomes 
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When  a smoke  layer  is  present,  the  plume  centerline  scalar  concentration  must  be  corrected  for 
the  entrained  scalar  as  the  smoke  plume  passes  through  the  smoke  layer.  This  can  be 
approximated  by  replacing  the  fire  with  a substitute  source  and  requiring  that  the  plume 
produced  by  the  substitute  source  is  entirely  in  the  smoke  layer.  The  location  and  strength  of 
the  substitute  source  may  be  determined  by  requiring  that  the  scalar  mass  flux  at  the  layer 
interface  be  conserved  and  that  the  plume  centerline  scalar  concentration  be  continuous  across 
the  layer  interface.  The  requirement  that  the  plume  centerline  scalar  concentration  be 
continuous  across  the  layer  interface  is  expressed  by 


16 


where  Cl  is  the  scalar  concentration  in  the  layer  and  the  subscripts  1 and  2 refer  to  the  plume 
centerline  scalar  concentration  at  the  layer  interface  in  the  lower  layer  and  the  plume  centerline 
scalar  concentration  in  the  smoke  layer  at  the  equivalent  layer  interface  height.  Using  equation 
15,  the  substitute  source  is 
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The  conservation  of  scalar  mass  flux  at  the  boundary  may  be  expressed  by 
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Substituting  equations  2 and  3 into  18  and  integrating  yields  with  the  substitution  of  equations 
5,  6,  and  14  the  location  of  the  substitute  source  with  respect  to  the  layer  interface. 
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The  requirement  that  the  distance  to  the  ceiling  from  the  interface  must  be  identical  for  both 
calculations  is  satisfied  by 

H2~Z2  = H]  — Z, 
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The  plume  centerline  scalar  concentration  at  the  ceiling  can  then  be  calculated  using  the 
substitute  source  by 
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The  Turning  Region 

The  next  step  is  to  obtain  the  smoke  or  gas  concentration  in  the  ceiling  jet  using  the  calculated 
scalar  concentration  in  the  plume.  Following  Alpert ' s derivation6  and  equating  the  mass  flux  in 

*7 

the  plume  to  mass  flux  at  the  start  of  the  ceiling  jet,  Yamauchi'  developed  an  equation  which 
related  the  maximum  smoke  concentration  in  the  plume  at  the  ceiling,  Cspo  to  the  average  smoke 
concentration  at  the  start  of  the  ceiling  jet,  Cs>aVe  in  terms  of  the  Gaussian  width  ratio  A for  the 
velocity  and  temperature  profiles  in  the  plume.  Assuming  that  the  same  ratios  work  for  gas  as 
well  as  smoke,  the  scalar  concentration  for  gas  or  smoke  is  given  by 

A2 

CSMVe  = Yri{C^~CL)+CL  22 

Assuming  that  the  scalar  concentration  of  gas  or  smoke  in  the  unconfined  ceiling  jet  can  be 
represented  by  a half  Gaussian  profile,  the  maximum  scalar  concentration  in  the  ceiling  jet,  Cs  is 
given  by 


A2 

Cjr=  0.18H)  = V2  j—r(Csp0(H)- q)+  Cl  2 

where  A2  = 1.157,  r is  the  radial  distance  from  the  plume  centerline  and  H is  the  distance  from 
the  surface  of  the  fire  to  the  ceiling. 
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The  Ceiling  Jet  Region 

The  scalar  concentration  of  gas  or  smoke  in  the  ceiling  jet  may  be  calculated  from  the  scalar 
mass  flux  equation  by  integrating  over  the  vertical  dimension,  y. 

1 i 00 

--7- J[(ci(/-,y)-  CL)v(r,y)rdy\=  0 24 

? ar  0 

It  has  been  assumed  that  the  density  of  the  scalar  concentration  in  the  layer  is  constant, 
independent  of  position.  The  resulting  spatial  averages  yield  the  average  scalar  concentration  in 
the  ceiling  jet  as  a function  of  r as 


c»=(cTJ-cJ^+c, 


rvh 
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where  the  subscript  “e”  represents  the  location  where  the  ceiling  jet  forms (re- 0.1 8H),  h is  the 
average  thickness  of  the  ceiling  jet,  v is  the  average  ceiling  jet  velocity  and  r is  the  radial 
distance  from  plume  center.  Using  Alpert 7 s correlation  for  the  maximum  ceiling  jet  velocity 
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and  fitting  Alpert 7 s calculation6  for  ceiling  jet  thickness  to  a power  law  {h/H  ~ (r/H)0'4, 
0.18H<r<2.0H),  the  maximum  scalar  concentration  in  the  ceiling  jet  is  given  by 


QoW  “ (Qo(0  Cl) 
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where  the  ceiling  jet  is  assumed  to  be  a half  Gaussian.  Replacing  Cso(re)  using  Equations  21 
and  23,  the  maximum  scalar  concentration  at  a radial  location  from  plume  center  in  the  ceiling 
jet  for  r > 0.18  H is  given  by 
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With  the  appropriate  product  yield  fraction  in  equations  8 and  10,  equation  28  predicts  the 
concentration  of  smoke,  CO,  or  other  combustion  product  in  the  ceiling  jet. 

Experiment 

Experiments  were  conducted  to  test  the  capability  of  the  algorithm  to  predict  CO  and  smoke 
concentration  in  the  ceiling  jet.  The  experiments  were  conducted  in  a room  with  floor 
dimensions  of  3.15  m by  3.02  m.  A sand  burner  was  centered  in  the  room  with  the  height 
between  the  top  of  the  burner  and  the  ceiling  being  either  2.19  m,  1.50  m or  0.77  m.  Two 
burner  geometries  were  used.  The  first  burner  was  round  with  a diameter  of  0.085  m while  the 
second  burner  was  square  with  an  effective  diameter  of  0.194  m.  The  burners  were  operated 
with  heat  release  rates  sufficient  that  the  flame  heights  should  scale  as  the  2/5  power  of  the  heat 
release  rate  (HRR).11  The  flat,  smooth  ceiling  was  made  up  of  acoustic  ceiling  tile.  The  walls 
of  the  room  were  constructed  of  glazed  cinderblocks  except  for  one  that  was  constructed  of 
ceiling  tile. 

Instrumentation  included  thermocouples,  CO  detectors,  a laser  extinction  meter  to  measure 
obscuration  and  a tapered-element  oscillating  microbalance  (TEOM)  to  measure  smoke  mass 
concentration.  Temperature  measurements  in  the  ceiling  jet  were  conducted  using  bare  bead 
0.076  cm  thermocouples  positioned  1 cm  below  the  ceiling  at  radial  positions  shown  in  figure  1. 
A thermocouple  tree  located  60  cm  from  plume  center  with  thermocouples  located  vertically 
below  the  ceiling  at  1.0,  3.0,  6.0,  12.0,  24.0,  36.0,  50.0,  75.0  and  100.0  cm  were  used  to 
measure  the  depth  of  the  ceiling  jet  and  the  temperature  in  the  smoke  layer.  An  aspirated 
thermocouple  was  used  to  check  the  impact  of  radiation  on  the  bare  bead  thermocouples. 


The  CO  measurements  were  conducted  using  electrochemical  cells  removed  from  residential 
CO  detectors.  The  carbon  monoxide  sensors  inside  the  detectors  consisted  of  electrochemical 
cells  with  terminals  connected  to  a resistor.  The  voltage  drop  across  the  resistor  increased 
linearly  when  CO  was  detected.  The  electrochemical  cells  were  calibrated  by  exposing  them  to 
a controlled  amount  of  CO  gas  inside  the  fire  emulator/detector  evaluator  tunnel  (FE/DE)12. 

The  detectors  were  placed  in  the  empty  evaluation  tunnel  and  the  voltage  reading  corresponding 
to  0 volume  fraction  of  CO  was  recorded  for  120  seconds.  Then,  a mixture  of  CO  gas  and  air 
was  allowed  to  flow  into  the  tunnel.  The  mixture  had  a constant  volume  fraction  equal  to  30  x 
10'6  of  CO,  which  was  verified  by  a (Siemens  Ultramat*)  gas  analyzer  using  IR  absorption  to 
measure  the  gas  concentration.  The  CO  mixture  flowed  until  the  electrochemical  cells  reached 
a steady  voltage,  at  which  point  the  CO  mixture  continued  to  flow  while  60  seconds  of  voltage 
data  was  recorded.  The  CO  flow  was  then  turned  off  and  the  CO  cleared  from  the  tunnel  to 
confirm  that  the  cells  would  return  to  their  original  voltages.  A mean  value  for  the  voltage  at  a 
volume  fraction  of  CO  of  0 and  30  x 10'6  was  calculated  for  each  cell.  These  two  points  were 
used  to  scale  the  CO  from  the  voltage,  assuming  a linear  relationship,  which  had  been 
determined  in  earlier  testing.  The  housing  surrounding  the  CO  detector  would  accept  flow 
between  the  ceiling  and  3.2  cm  beneath  the  ceiling.  There  were  also  small  vents  located  at  the 
bottom  of  the  housing  about  8.9  cm  beneath  the  ceiling.  Figure  1 provides  the  layout  for  the 


* The  identification  of  any  commercial  product  or  trade  name  does  not  imply  endorsement  or  recommendation  by 
the  National  Institute  of  Standards  and  Technology. 
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detectors  with  detectors  one  through  six  mounted  on  the  ceiling  and  detector  seven  mounted  at 
1/3  of  the  distance  between  the  burner  and  the  ceiling.  Detectors  3 and  5 were  located  along 
the  same  radius  and  were  used  to  check  the  radial  dependence  of  the  CO  volume  fraction. 

Smoke  obscuration  measurements  were  conducted  using  a helium-neon  laser  mounted  such  that 
the  exposed  laser  beam  was  located  1.0  m from  the  plume  centerline,  6.35  cm  beneath  the 
ceiling  and  had  a path  length  of  0.5  m as  shown  in  figure  1.  By  measuring  the  decrease  in 
intensity  of  the  beam,  the  smoke  density  is  deduced  from  the  light  extinction  coefficient  using 

9 io  ^ ^ 

the  extinction  coefficient  per  unit  mass,  8.71  irT/g  . 

A tapered  element  oscillating  microbalance  (TEOM)  was  used  as  a second  method  of  measuring 
smoke  concentration.  The  sampling  tube  for  the  TEOM  was  positioned  4.4  cm  beneath  the 
ceiling  and  1.0  m from  the  plume  centerline. 

Two  fuels,  propane  (CsHg)  and  propene  (C3H6),  were  used  for  the  experiments.  Values  for  the 
heat  of  combustion,  radiation  fraction,  smoke  yield  and  CO  yield  for  these  fuels  are  taken  from 
the  literature  and  tabulated  below9 13 14 15 16.  The  calculations  were  done  using  the  values  of 
Tewarson9  or  using  values  obtained  from  the  NIST  cone  calorimeter.  Differences  in  the  yield 
fractions  and  heat  of  combustion  enter  the  calculations  linearly.  Differences  in  radiation 
fraction  will  have  a small  and  nonlinear  impact  on  the  calculations  as  the  depth  of  the  smoke 
layer  will  be  impacted  as  well  as  values  in  the  correlation  developed  above. 

Table  1.  Yield  fractions  for  smoke  and  CO,  radiation  fraction,  and  heat  of  combustion  for 
propane  and  propene. 


Ref# 

Ys 

y 

1 CO 

Radiation  Fraet. 

Heat  of  Comb. 

Propane/Propene 

Propane/Propene 

Propane/Propene 

Propane/Propene 

(g/g) 

( g/g) 

kJ/g 

9 

0.024/0.095 

0.005/0.017 

0.27/0.32 

43.7/40.5 

13 

0.021/0.074 

14 

0.010/0.069 

0.28/0.40 

15 

0.27/0.38 

16 

.0045/x 

0.27/x 
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Analysis 


The  predictive  capability  of  the  ceiling  jet  algorithm  is  compared  with  measurements  of  CO  and 
smoke  in  the  ceiling  jet.  CO  measurements  were  also  conducted  in  the  smoke  layer.  The 
locations  of  the  ceiling  jet  and  layer  interface  were  determined  using  the  temperature 
measurements  from  the  thermocouple  tree  located  60  cm  from  the  fire  center.  The  ceiling  jet 
can  be  identified  as  the  region  near  the  ceiling  in  which  the  temperature  is  about  uniform.  For 
the  experiments  described  in  this  paper,  it  can  be  seen  from  figs.  2 through  6 that  the  ceiling  jet 
boundary  develops  rapidly  and  extends  down  to  between  6 and  12  cm  beneath  the  ceiling  in  all 
cases. 

CO  Analysis 

The  CO  measurements  in  the  ceiling  jet  and  smoke  layer  were  made  using  commercial  CO 
detectors.  The  response  characteristics  of  the  CO  detectors  must  be  known  in  order  to  compare 
the  detector  measurements  with  the  predictions  of  the  ceiling  jet  algorithm.  Since  the  fires  used 
in  the  experiment  were  small,  delay  times  must  also  be  included  to  account  for  the  time  for  the 
ceiling  jet  flow  to  reach  the  detector  and  the  time  required  to  form  a smoke  layer.  The  2.19  m 
propene  experiment  will  be  used  to  demonstrate  the  importance  of  these  factors. 

Figure  7 provides  a comparison  of  the  ceiling  jet  algorithm  predictions  with  the  experimental 
data.  Transit  times  and  detector  characteristics  are  not  included  in  the  modeling.  The  zone 
model  calculations  assume  that  the  ceiling  jet  is  fully  formed  and  that  the  layer  begins  to 
develop  at  the  start  of  the  experiment.  The  zone  model  JET  was  used  to  determine  the  layer 
height  and  temperature  for  the  calculations.  The  CO  volume  fraction  for  the  layer  was 
calculated  by  integrating  the  CO  production  over  time  and  dividing  by  the  layer  volume.  It  was 
assumed  that  there  was  no  loss  of  CO  from  the  layer.  As  seen  in  figure  7,  ignoring  the  transit 
times  as  well  as  the  detector  characteristics  produce  a substantial  over  prediction  of  the 
measured  CO  volume  fraction. 

The  detector  response  time  for  the  layer  and  layer  delay  time  can  be  included  in  the  calculation 
by  using  the  measured  CO  volume  fraction  for  the  layer  in  the  calculation  as  shown  in  figure  8. 
Since  the  zone  model  approximation  is  still  used  to  calculate  the  ceiling  jet  volume  fraction,  the 
calculated  value  for  the  ceiling  jet  starts  at  a nonzero  value  that  remains  constant  until  the 
smoke  layer  begins  to  form.  The  measured  value  for  the  ceiling  jet  begins  to  rise  at  30  s after 
the  start  of  the  experiment  and  approaches  the  calculated  value  for  the  ceiling  jet.  At  60  s after 
the  start  of  the  experiment,  the  CO  detector  in  the  layer  begins  to  respond  which  marks  the 
beginning  of  the  impact  of  the  smoke  layer  on  the  ceiling  jet.  Using  the  measured  value  of  the 
layer  CO  volume  fraction,  the  calculated  value  of  the  CO  volume  fraction  in  the  ceiling  jet  is 
nearly  identical  to  the  measured  value. 

In  order  to  model  the  response  of  the  CO  detectors  in  the  ceiling  jet  and  layer,  an  experimental 
effort  was  undertaken  to  develop  a detector  algorithm.  The  FE/DE  was  used  to  produce  a 
constant  velocity,  constant  CO  concentration  environment.  Smoldering  cotton  wicks  were  used 
to  supply  the  constant  source  of  CO.  The  detector  was  isolated  from  the  CO  environment  by 
placing  a tight  fitting  can  over  it.  Once  the  CO  concentration  near  the  detector  was  steady,  the 
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can  was  removed  and  the  detector  response  measured.  The  amount  of  time  required  for  the 
detector  to  begin  to  respond  once  it  was  exposed  to  CO  was  quite  short  and  taken  to  be  1.0  s. 
The  response  of  the  detector  could  be  simulated  using  the  fitting  function  y = l-(l  + r/  T)e~t,T . 
This  equation  is  consistent  with  a two-step  diffusion  process  with  the  time  constant  x being  the 
same  for  each  step.  The  gas  flow  velocity  was  varied  from  0.042  m/s  to  0.22  m/s  with  the 
detector’s  value  of  x varying  from  19.7  s to  27.9  s as  the  flow  velocity  decreased. 

The  ceiling  jet  and  detector  algorithms  for  smoke  and  gas  were  added  to  the  zone  model  JET 
and  the  model  predictions  were  tested  against  the  data  collected  with  the  CO  detectors.  The 
general  equation  set  used  to  represent  the  two-step  diffusion  process  for  the  detectors  is 


dXf  = [X,(,-Sr)-X(t)]/T 
dt 

29 

dX:{,)=[x(tyxs(t)Vr 

dt 

30 

St  - au~( 

31 

T = bu~d 

32 

X(t)  is  the  smoke  or  gas  mass  density  in  the  sensing  chamber,  Xe  is  the  smoke  or  gas  mass 
density  at  the  location  of  the  detector,  Xs  is  the  gas  mass  density  at  the  electrochemical  sensor, 
ue  is  the  flow  velocity  outside  of  the  detector  housing,  St  is  the  dwell  or  offset  time  for  the  flow 
to  travel  through  the  detector  to  the  sensing  chamber,  x is  the  mixing  time  in  the  sensing 
chamber,  and  a,  b,  c,  and  d are  determined  from  the  response  time  fits  of  the  detector.  For  flow 
velocities  above  0.1  m/s,  5x  and  x become  less  sensitive  to  velocity  variations  for  the  detectors 
tested  and  can  be  treated  as  constants,  which  was  the  case  for  these  tests. 

The  CO  mass  fraction  produced  by  the  small  burner  was  measured  using  the  NIST  cone 
calorimeter  and  found  to  be  0.0136.  That  is  slightly  smaller  than  the  value  0.017  given  in  table 
1.  The  value,  0.0136,  was  used  in  the  calculations. 

Two  propene  experiments  that  use  the  small  burner  with  a HRR  of  2.5  KW  and  have  ceiling 
heights  of  1.5  m and  2.19  m were  modeled  using  JET.  Due  to  the  changing  CO  concentrations 
in  the  ceiling  jet  and  layer,  the  CO  detectors  were  modeled  using  the  two  first  order  differential 
equations  that  simulated  the  fitting  function  for  constant  CO  exposure  found  for  the  detector. 
Time  constants  of  19.7  s and  27.9  s were  used  to  simulate  the  response  of  the  detector  in  the 
ceiling  jet  and  layer  respectively.  Transit  times  at  the  beginning  of  the  experiment  were 
measured  using  the  response  of  thermocouples  located  at  the  positions  of  the  CO  detectors  in 
the  ceiling  jet  and  layer.  For  the  2.19  m experiment,  the  transit  times  were  4 s to  the  ceiling  jet 
CO  detector  at  0.6  m and  28  s to  CO  detector  7 located  where  the  layer  should  form.  For  the  1.5 
m experiment,  the  measured  transit  times  were  3 s and  39  s to  these  detector  positions.  The  39  s 
time  was  difficult  to  estimate  as  the  temperature  increase  in  the  layer  was  quite  small  initially  at 
the  position  of  the  CO  detector  and  could  be  confused  with  the  slow  temperature  rise  of  the 
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thermocouples  responding  to  radiation  from  the  fire.  Transit  times  were  also  estimated  using 
the  velocity  equations  for  the  plume  and  ceiling  jet,  equations  6 and  26,  with  the  time  to  the 

ceiling  jet  detector  and  layer  detector  calculated  integrating  over  the  path  using  t = ^dxl  v 

where  v is  the  velocity.  It  was  assumed  that  the  return  flow  velocity  could  be  roughly 
approximated  using  equation  26  over  these  distances.  The  resulting  calculations  gave  values  of 
2 s and  39  s for  the  2.19  experiment  and  1 s and  45  s for  the  1.5  m experiment.  The  39  s and 
45  s times  are  the  times  for  the  ceiling  jet  to  flow  from  0.2  H to  the  wall  and  back  to  the  plume 
center  rather  than  to  the  0.6  m position. 

The  transit  times  to  the  ceiling  jet  detector  and  layer  detector  used  in  the  modeling  were  3 s and 
39  s for  the  1.5  m experiment  and  4 s and  28  s for  the  2.19  m experiment.  The  formation  of  the 
layer  was  delayed  in  the  model  using  the  layer  transit  time.  The  layer  begins  to  form  once  the 
ceiling  jet  strikes  the  walls  and  begins  to  turn  back  toward  the  plume.  The  ceiling  jet  conditions 
will  be  impacted  by  this  flow  once  the  flow  has  reached  the  vicinity  of  the  plume. 

Figures  9 and  10  provide  the  results  of  the  simulations  using  JET  and  are  compared  with  the 
measured  values  for  the  CO  volume  fraction  using  detectors  in  the  ceiling  jet  and  layer  at  0.6  m 
from  the  plume  centerline.  A small  baseline  displacement  of  roughly  -1.0  x 1CT6  in  the  volume 
fraction  occurred  for  the  ceiling  jet  detector  when  the  ceiling  jet  flow  interacted  with  the 
detector.  The  zero  level  for  both  figures  was  taken  as  the  bottom  of  the  dip.  If  the  baseline 
returned  to  its  normal  location  after  the  initial  thermal  insult,  the  volume  fraction  shown  on  the 
curves  for  the  ceiling  jet  measurement  should  be  decreased  by  1.0  x 1CT6.  No  baseline  shift  was 
observed  for  the  detectors  in  the  smoke  layer. 

The  model  predictions  follow  the  increasing  trend  of  the  measured  volume  fraction  for  the  first 
80  s of  the  experiment.  At  later  times  the  measured  values  rise  more  slowly  than  the  predicted 
values.  The  predicted  and  measured  values  of  CO  volume  fraction  agree  to  within 
approximately  25  % or  better  over  the  experimental  range.  The  fall-off  in  the  agreement 
between  the  model  and  experiments  at  the  later  times  could  result  from  the  slow  venting  of  the 
smoke  layer  as  the  model  did  not  include  losses  in  the  calculations. 

Comparisons  of  the  measured  values  of  CO  using  the  large  burner  with  a fire  size  of  7.6  kW  and 
ceiling  heights  of  1.5  m and  2.19m  are  given  in  figures  1 1 and  12.  The  CO  mass  yield 
measured  for  the  small  propene  burner  was  used  to  represent  the  large  propene  burner,  as  the 
large  burner  was  too  big  to  be  tested  in  the  cone  calorimeter.  The  baseline  displacement  for  the 
ceiling  jet  detectors  were  larger  for  these  experiments  with  two  curves  given  for  the  ceiling  jet 
volume  fraction,  one  that  uses  the  baseline  prior  to  the  start  of  the  experiment  and  the  other  that 
takes  the  bottom  of  the  baseline  displacement  as  the  zero  point.  The  detectors  in  the  smoke 
layer  did  not  show  baseline  displacements.  The  predicted  trend  of  the  CO  volume  fraction  in 
the  smoke  layer  follows  the  measurements  for  both  heights.  The  predicted  trend  for  the  ceiling 
jet  follows  the  measurements  for  the  2.19  m height  but  increases  with  a smaller  slope  for  the  1.5 
m height.  The  measured  values  for  CO  in  the  smoke  layer  agreed  with  the  prediction  but  in  the 
ceiling  jet,  the  predicted  values  lay  below  the  measured  values  with  the  magnitude  of  the 
discrepancy  dependent  on  the  amount  of  baseline  shift. 
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Comparisons  using  propane  with  the  small  burner  at  a fire  size  of  2.7  kW  are  given  in  figure  13. 
There  is  a baseline  displacement  of  2 x 10'6  in  volume  fraction  for  the  detector  in  the  ceiling  jet 
and  the  comparisons  given  in  the  figure  include  the  2 x 10"6  baseline  shift.  The  predicted  trend 
for  the  CO  in  the  ceiling  jet  follows  the  experimental  trend  although  in  the  layer,  the  measured 
trend  does  not  follow  the  predicted  trend.  It  is  possible  that  the  CO  volume  fraction  in  the  layer 
is  so  small  that  the  measurement  is  more  noise  than  signal. 

Smoke  Analysis 

Figures  14-16  display  the  comparisons  between  the  measured  smoke  density  and  the  calculated 
smoke  density  for  the  upper  layer  and  the  maximum  value  in  the  ceiling  jet  using  JET.  For  the 

2.5  kW  fire,  the  predictions  either  lie  below  or  between  the  predictions  for  the  density  in  the 
layer  and  ceiling  jet  with  the  2.19  m measurement  showing  a sudden  increase  in  measured 
density  at  about  170  s into  the  experiment.  This  sudden  increase  in  density  was  observed  in  the 

7.6  kW  experiment  and  could  be  a temperature  problem  with  the  diode  detector  or  laser.  The 
slope  of  the  measured  values  did  not  agree  with  the  calculated  values  for  the  7.6  kW  experiment 
but  were  much  better  for  the  2.5  kW  experiments.  The  measured  values  from  the  TEOM  are 
consistently  below  the  predicted  values  as  demonstrated  in  figure  16.  Measurements  were  not 
available  for  the  smoke  density  in  the  layer. 

Radial  Dependence  Analysis 

The  radial  dependence  of  the  ceiling  jet  algorithm  was  tested  using  CO  detectors  3 and  5. 
Detector  5 was  located  60  cm  from  plume  center  and  detector  3 was  located  150  cm  from  plume 
center  along  the  same  radius.  Based  on  the  distance  difference  of  the  two  detectors,  detector  3 
should  give  a signal  0.59  smaller  than  detector  5 if  the  radial  decrease  of  the  CO  concentration 
in  the  ceiling  jet  follows  equation  28.  Figures  17  and  18  compare  the  signal  recorded  by 
detector  3 with  the  signal  from  detector  5 multiplied  by  0.59  for  the  2.5  kW  propene  fires  at 
2.19  m and  1.5  m ceiling  heights.  The  curves  agree  within  the  measurement  uncertainty  of  the 
experiment  that  is  estimated  to  be  1 x 10~6  based  on  the  uncertainty  of  the  zero  measurement  for 
the  two  detectors.  No  correction  was  made  for  time  of  flight  between  the  two  detectors  but  if  a 
correction  were  made  the  agreement  would  improve,  as  the  curve  for  detector  3 would  move  to 
the  left  to  compensate  for  the  increased  time  for  the  hot  gas  to  reach  that  detector. 

Conclusion 

Algorithms  to  predict  gas  and  smoke  concentrations  in  the  plume  and  ceiling  jet  have  been 
developed  and  tested  against  experiments.  The  zone  model  JET  was  modified  to  include  the 
algorithms  and  was  used  to  generate  the  model  predictions.  The  comparison  of  prediction 
versus  measurement  for  CO  in  the  ceiling  jet  demonstrated  the  importance  of  including  the 
detector  response  and  the  development  time  for  the  smoke  layer  in  the  calculations. 

The  CO  sensing  elements  in  the  ceiling  jet  responded  to  the  changing  temperature  of  their 
environment  by  shifting  the  zero  baseline  at  the  start  of  the  experiment.  This  shift  in  zero 
baseline  introduced  an  uncertainty  in  the  measured  CO  volume  fraction  that  was  more 
pronounced  for  propane  with  its  low  level  of  CO  production  and  for  the  7.6  kW  propene 
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experiments  with  their  higher  temperatures.  It  would  be  desirable  to  measure  the  CO  volume 
fraction  by  an  independent  method  to  check  the  consistency  of  this  measurement. 

The  experiments  used  for  the  comparison  of  the  algorithm  predictions  with  the  smoke 
concentration  measurements  bear  repeating.  The  discrepancy  in  the  smoke  concentration 
measurements  using  the  TEOM  and  laser  attenuation  also  needs  to  be  investigated.  The  smoke 
concentration  values  provided  by  the  TEOM  were  lower  than  expected.  The  laser  signals  were 
noisy  which  made  it  difficult  to  determine  whether  the  laser  optics  were  being  coated  by  smoke 
during  the  experiment.  The  impact  of  temperature  effecting  the  laser  measurements  also  needs 
to  be  investigated. 

Additional  measurements  of  the  smoke,  CO,  and  temperature  as  a function  of  depth  need  to  be 
done  in  order  to  verify  that  the  structure  of  the  ceiling  jet  is  similar  to  the  structure  determined 
by  temperature  measurements. 

In  summary,  ceiling  jet  algorithms  have  been  developed  for  predicting  smoke  and  CO 
concentrations  in  the  presence  of  a growing  layer.  These  algorithms  can  be  used  with 
appropriate  detector  algorithms  to  predict  the  activation  of  smoke  and  CO  detectors  located  in 
the  ceiling  jet.  The  new  algorithms  have  been  incorporated  in  the  zone  model  JET  and  will 
become  available  when  the  next  version  of  JET  is  released. 
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Figure  1 Detector  locations  for  the  experiments.  The  solid  circles  are  thermocouples,  the  rings 
are  CO  detectors,  and  the  location  of  the  laser  beam  is  shown  on  the  right.  The  TEOM  inlet  was 
located  just  below  the  laser  beam  at  100  cm  from  the  burner.  The  CO  detector  tree  consisted  of 
a detector  in  the  ceiling  jet  with  additional  detectors  located  at  1/3  and  2/3  the  distance  between 
the  burner  and  the  ceiling.  The  thermocouple  tree  had  thermocouples  located  at  (1,  3,  6,  12,  24, 
36,  50,  75,  and  100)  cm  below  the  ceiling. 
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Figure  2 Excess  temperature  measurements  for  propene  at  60  cm  from  the  axis  of  the  small 
burner  (flow  rate  of  2.0  L/min)  that  was  located  1.5  m below  the  ceiling. 
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Figure  3 Excess  temperature  measurements  for  propene  at  60  cm  from  the  axis  of  the  large 
burner  (flow  rate  of  6.0  L/min)  that  was  located  1.5  m below  the  ceiling. 
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Figure  4 Excess  temperature  measurements  for  propene  at  60  cm  from  the  axis  of  the  small 
burner  (flow  rate  of  2.0  L/min)  that  was  located  2.19  m below  the  ceiling. 
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Figure  5 Excess  temperature  measurements  for  propane  at  60  cm  from  the  axis  of  the  small 
burner  (flow  rate  of  1.88  L/min)  that  was  located  2.19  m below  the  ceiling. 
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Figure  6 Excess  temperature  measurements  for  propane  at  60  cm  from  the  axis  of  the  small 
burner  (flow  rate  of  1.88  L/min)  that  was  located  1.5  m below  the  ceiling. 
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Figure  7 Comparison  of  measured  and  calculated  CO  volume  fraction  for  the  2.5  kW  small 
burner  experiment  at  a height  of  2.19  m.  Delay  times  and  detector  response  functions  are  not 
included  in  the  calculation. 
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Figure  8 CO  volume  fraction  measured  in  the  ceiling  jet  and  layer  (mea)  and  calculated  for  the 
ceiling  jet  (cal).  The  fire  was  propene  at  2.5  kW  (small  burner)  with  the  top  of  the  burner 
located  2.19  m beneath  the  ceiling.  The  measured  value  of  the  CO  volume  fraction  was  used  in 
the  calculation  of  the  ceiling  jet  molar  fraction. 
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Figure  9 CO  volume  fraction  measured  in  the  ceiling  jet  and  layer  (mea)  and  calculated  for  the 
ceiling  jet  (cal).  The  fire  was  propene  at  2.5  kW  (small  burner)  with  the  top  of  the  burner 
located  2.19  m beneath  the  ceiling.  Ceiling  jet  transit  time,  layer  formation  delay  time  and 
detector  response  characteristics  are  included  in  the  calculations. 
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Figure  10  CO  volume  fraction  measured  in  the  ceiling  jet  and  layer  (mea)  and  calculated  for  the 
ceiling  jet  (cal).  The  fire  was  propene  at  2.5  kW  (small  burner)  with  the  top  of  the  burner 
located  1.5  m beneath  the  ceiling.  Ceiling  jet  transit  time,  layer  formation  delay  time  and 
detector  response  characteristics  are  included  in  the  calculations. 
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Figure  11  CO  volume  fraction  measured  in  the  ceiling  jet  and  layer  (m  and  mea)  and  calculated 
for  the  ceiling  jet  and  layer  (cal).  The  ceiling  jet  measurements  are  plotted  using  the  zero 
baseline  as  the  bottom  of  the  signal  dip  (dp)  and  as  the  baseline  at  the  start  of  the  experiment 
(ndp).  The  fire  was  propene  at  7.6  kW  (large  burner)  with  the  top  of  the  burner  located  2.19  m 
beneath  the  ceiling.  Ceiling  jet  transit  time,  layer  formation  delay  time  and  detector  response 
characteristics  are  included  in  the  calculations. 
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Figure  12  CO  volume  fraction  measured  in  the  ceiling  jet  and  layer  (m  and  mea)  and  calculated 
for  the  ceiling  jet  (cal).  The  ceiling  jet  measurements  are  plotted  using  the  zero  baseline  as  the 
bottom  of  the  signal  dip  (dp)  and  as  the  baseline  at  the  start  of  the  experiment  (ndp).  The  fire 
was  propene  at  7.6  kW  (large  burner)  with  the  top  of  the  burner  located  1.5  m beneath  the 
ceiling.  Ceiling  jet  transit  time,  layer  formation  delay  time  and  detector  response  characteristics 
are  included  in  the  calculations. 
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Figure  13  CO  volume  fraction  measured  in  the  ceiling  jet  and  layer  (mea)  and  calculated  for  the 
ceiling  jet  (cal).  The  fire  was  propane  at  2.7  kW  (small  burner)  with  the  top  of  the  burner 
located  2.19  m beneath  the  ceiling.  Ceiling  jet  transit  time,  layer  formation  delay  time  and 
detector  response  characteristics  are  included  in  the  calculations. 
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Figure  14  Smoke  density  measured  in  the  ceiling  jet  (laser  m)  and  calculated  for  the  ceiling  jet 
(CJ  cal)  and  smoke  layer  (Lyr  cal).  The  fire  was  propene  at  2.5  kW  (small  burner)  with  the  top 
of  the  burner  located  1.5  m beneath  the  ceiling.  Ceiling  jet  transit  time  and  layer  formation 
delay  time  are  included  in  the  calculations. 
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Figure  15  Smoke  density  measured  in  the  ceiling  jet  (laser  and  TEOM)  and  calculated  for  the 
ceiling  jet  (CJ  cal)  and  smoke  layer  (Lyr  cal).  The  fire  was  propene  at  2.5  kW  (small  burner) 
with  the  top  of  the  burner  located  2.19  m beneath  the  ceiling.  Ceiling  jet  transit  time  and  layer 
formation  delay  time  are  included  in  the  calculations. 
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Figure  16  Smoke  density  measured  in  the  ceiling  jet  (laser  m)  and  calculated  for  the  ceiling  jet 
(Cj  cal)  and  smoke  layer  (Lyr  cal).  The  fire  was  propene  at  7.6  kW  (small  burner)  with  the  top 
of  the  burner  located  1.5  m beneath  the  ceiling.  Ceiling  jet  transit  time  and  layer  formation 
delay  time  are  included  in  the  calculations 
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Figure  17  Radial  dependence  comparison  using  CO  detectors  5 and  3 located  in  the  ceiling  jet 
at  60  cm  and  150  cm  respectively.  Detector  5 has  had  its  zero  corrected  by  0.6  and  is  scaled  to 
0.59  which  is  the  predicted  scaling  by  the  ceiling  jet  algorithm.  A 2.5  kW  propene  fire  was 
used  (small  burner)  at  a distance  2.19  m below  the  ceiling. 
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Figure  18  Radial  dependence  comparison  using  CO  detectors  5 and  3 located  in  the  ceiling  jet 
at  60  cm  and  150  cm  respectively.  Detector  5 has  had  its  zero  corrected  by  0.6  and  is  scaled  to 
0.59  which  is  the  predicted  scaling  by  the  ceiling  jet  algorithm.  A 2.5  kW  propene  fire  was  used 
(small  burner)  at  a distance  1.5  m below  the  ceiling. 
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